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ABSTRACT 


A theoretical analysis ami an experimental investigation were conducted to assess 
the feasibility of developing a spinning wave heat engine. Such an engine would utilize a 
large amplitude traveling acoustic wave rotating around a cylindrical chamber, and it 
should not suffer from the inefficiency, noise, and intermittent thrust which characterizes 
pulse jet engines. The objective of this investigation was to determine whether an 
artificialiy driven large amplitude spinning transverse wave could induce a steady flow of 
air through the combustion chamber under cold flow conditions. In the theoretical analysis 
the Maslen and Moore perturbation technique was extended to study flat cylinders 
(pancake geometry) with completely open side walls and a central opening. In the parallel 
experimental study, a test model was used to determine resonant frequencies and radial 
pressure distributions, as well as oscillatory and steady flow velocities at the inner and 
outer peripheries. The experimental frequency was nearly the same as the theoretical 
acoustic value for a model of the same outer diameter but without a central hole. 
Although the theoretical analysis did not predict a steady velocity component, 
simultaneous measurements of hotwire and microphone responses have shown that the 
spinning wave pumps a mean flow radially outward through the cavity. This conclusion has 
been verified by flow visualization using both a sensitive tufted wand and smoke as flow 
tracers. These results indicate that further development and testing of the spinning wave 
engine concept is warranted. 
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INTRODUCTION 


This report is a summary of work conducted under NASA Grant No. 
NAG 3-96 during the period September 8» 1980 through September 7t 1981. 
This project is entitled "Development of a Spinning Wave Heat Engine" and 
the NASA Technical Officer is B. R. Phillips at the Lewis Re iearch Center. 

The purpose of this research program is to provide basic information 
needed for the development of a spinning wave heat engine originally 
conceived by M. F. Heidmann at the NASA Lewis Research Center. The 
spinning wave heat engine concept is based on the presence of a large 
amplitude traveling acoustic wave rotating around a cylindrical chamber* 
The spinning wave engine is conceptually similar to the pulse jet engine, but 
fundamental differences between the two engines suggest that the spinning 
wave concept can eliminate or alleviate the problems of inefficiency, noise, 
and intermittent thrust that caused the demise of pulse jet engine 
development. On the other hand, the spinning wave engine retains the 
characteristics of simplicity, low cost, low weight and absence of rotating 
parts that stimulated the development of pulse jets in the early 1950's. 

As a propulsive unit by itself, a spinning wave engine is not a 
competitor of current engines for commercial CTOL and V/STOL aircraft. 
Effective pressure ratio is relatively low (less than 3:1) and, therefore, 
specific fuel consumption for such applications would be relatively high. The 
concept, however, remains attractive for many applications such as military, 
powered glides, auxiliary power, precompression stage for other engine 
cycles, and ground-effect vehicles. Within its pressure ratio range, the cycle 
efficiency of a spinning wave engine can be higher than that of other cycles 


because heat is essentially added at con$tant volume. 

In order for the spinning wave engine to function, the spinning 
acoustic wave must induce a steady flow of air through the combustion 
chamber. Such a steady flow is a property of nonlinear waves (i.e., large 
amplitude) in which the oscillating flow in one direction is greater than that 
in the opposite direction and a net through flow result$. In the nonlinear 
analysis by Maslen and Moore\ large amplitude standing and traveling 
acoustic waves in a closed circular cylinder were studied using a 
perturbation technique. Maslen and Moore's results indicate that a small 
depression of mean pressure occurs at the center and a small elevation of 
mean pressure occurs at the periphery when a large amplitude traveling first 
tangential (IT) wave is present. However, no net through flow can occur for 
this geometry because there are no inlet and outlet openings. When viscosity 
is included in the analysis, a steady wheel flow is predicted which rotates in 
a direction opposite to the direction of wave travel for low order modes. 
This wheel flow, however, does not contribute to a net flow through the 
chamber. 

In order to obtain a steady flow through the combustor, an inlet and 
outlet must be provided. The results of Maslen and Moore for the spinning IT 
mode appear to imply that if openings are provided at the center and at the 
periphery of the cylinder, a net inflow will occur at the center and an 
outflow will occur at the periphery. Providing the openings, however, will 
modify the boundary conditions at the center and the periphery so that 
Maslen and Moore's analysis no longer applies. In the most extreme case the 
side wall of the cylinder is completeiy open and the boundary condition' at 


the wall requires (approximately) that the pressure oscillations vanish there 
(in contrast to the closed cylinder in which the radial velocity component 
vanishes at the wall). Similarly^ pr^wlding a central inlet Is expected to 
modify the structure of the spinning transverse modes. 

In this investigation tho acoustic generation cf a steady through flow 
in the absence of combustion is determined by extending the Masien and 
Moore analysis to include the effects of inlet and outlet openings. This work 
is divided into two parts. In the first part the Masien and Moore technique is 
used to study flat cylinders (pancake geometry) with completely open side 
walls and a central opening. In the second part an experimental apparatus is 
used to verify the analytical results. The remainder of this report 
summarizes the analytical and experimental work done under this project. 
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THEORETICAL ANALYSIS 

The objective of the theoretical analysis described herein is to 
determine whether a large amplitude spinning acoustic wave can pump a 
steady flow through a flat cylindrical chamber which is open on the 
periphery and at the center. Such an annular pancake cavity is formed by 
two parallel, closely spaced annular disks as shown in Figure 1, To minimize 
losses due to acoustic radiation and to ensure a two-dimensional acoustic 
solution, the spacing h between the disks must be small compared to the 
outer radius b; that is h/b <<1. Under these conditions the pressure at the 
open boundaries is approximately equal to the ambient pressurej that is, as a 

first approximation the pressure perturbation must vanish at the outer edge 

^ # 

(r s b) and inner edge (r = a) of the annular plates. This is analogous to the 
open-end boundary condition used in the analysis of organ pipes. In contrast, 
the nonlinear analysis of Maslen and Moore^^^ was applied to hard walled 
cylinders for wnich the radial velocity component must vanish at the outer 
boundary (there is no inner boundary in this case). In fact for the first 
tangential spinning mode, the maximum pressure fluctuation occurs at the 
hard walled boundary in the Maslen and Moore analysis. 

Although the boundary conditions are different, the governing 
differential equations are the same for both geometries. Therefore the 
approach taken in this analysis is to follow closely the Maslen and Moore 
procedure, substituting the appropriate boundary conditions where 
necessary. This leads to some serious difficulties, which have not been 
satisfactorily resolved. Nevertheless it is believed that the solutions 
obtained in this manner are instructive and shed some light on the expected 



f t 


6 


behavior of the spinning wave engine. The exploration of alternative 
methods of analysis^ such as the Method of Weighted Residuals, was 
beyond the scope of this one-year project. 

Modification of Maslen and Moore Analysis 

Basic Procedure. In the Maslen and Moore analysis it is assumed that 
the flow is isentropic and irrotational so that a velocity potential, exists 
such that ^ s u. The energy and state equations yield the isentropic flow 
relation, P s . The continuity and momentum equat'ons yield the 
following equatiof/ ;ior the velocity potential} 

«P • V V) V + J VV • V(7<p • V«p) 

+ (Y-1) Vt 7^<p (1) 

where y is the specific heat ratio and the pressure is related to <p by: 

^ r 

1 - p Y , (y-l) 2 J 

The above equations involve only dimensionless quantities defined as follows: 


p 


nondimensional pressure, P /P^ 

t 


nondimensional time, cp t 

u,v,w 

s 

nondimensional axial, radial, and tangential velocities 

# # # 

(cylindrical coordinates), u /c^, v w /c^ 

U 

= 

nondimensional velocity vector;, 

X 


nondimensional axial coordinate, mx /c 

o 

O’ 


# 

nondimensional radial coordinate, wr /c^ 
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9 » itr^gle in cylindrical coordinates 

# ' 

p * nondimensional densityr p /p ^ 

where u is the angular frequency (dimentional)i c^ is the velocity of sound in 
the absence of the wave> the subscript o refers to dimensional flow 
properties in the absence of the wave, and the superscript* refers to 
dimensional properties. 

In the perturbation procedure used by Maslen and Moore the velocity 
potential and pressure are expanded in terms of an amplitude parameter c , 
thus 




(3) 

and 


(4) 


Substituting Eqs. (3) and (4) into Eqs. (1) and (2)» collecting and grouping 
terms according to powers of e , and equating the resulting coefficients of 
each power of e to zero gives equations to be solved successively for the 
first) second and third order potentials and pressures. 

m 

First Order Solutions . Equating the coefficient of « in Eq. (1) to 
zero gives the wave equation for the first order potential 

, ( 5 ) 

while Eq. (2) yields: 


( 6 ) 


• t * 
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Spinning wave solutions of Eq. (5) were obtained by the method of 
separation of variables and were required to satisfy the boundary conditions 
P' ^ X 0 at r s a (inner) and r s b (outer) boundaries. These acoustic 
solutions for the velocity potential have the formt 



V^( Pni^r)! sin (t nS ) (7) 


where 3^ and are Bessel functions of the first and second kinds, 
n n ' 

respectively, of order n. The parameter v is the ratio of inner radius to 
outer radius, v = a/b, and the eigenvalue 8 is the k' solution Of the 
equation 


3_(x)Y ( V x) - Y (x)3 ( V x) B 0. (8) 

n n n n 

The velocity potential is expressed in terms of the dimensionless radius 
# 

r s r /b, the angular coordinate 6 , and the dimensionless time t. The actual 
frequency f is then related to the eigenvalue by 


The acoustic pressure is then given by 
P^*^ B - Y R„( o) cos (t + n0 ) 


( 10 ) 


* * L 
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where R^(o^ ) is the radial acoustic eigenfunction given by 

3 J 

Rja) = Jj«) - JL-M. Yjoi) 

n n Y(vp .) " 

n nk 

and a = » r*/c^ = »„!< r*/'’ = S n^r- 
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( 11 ) 


Detailed calculations of the radial eigenfunctions, eigenvalues and 
pressure fields have been carried out for the principal mode of interest, the 
spinning first tangential mode (n ? 1, k = 1). A plot of the eigenvalue ^ as 
a function of v is shown in Figure 2, which implies that the frequency 
increases as the diameter of the central opening increases. For v< 0.1 
Figure 2 also shows that the frequency is close to the value ( Pjj = 3.83171) 
obtained when no central opening is present (i.e,, two closely spaced disks). 
The radial dependence of is shown in Figure 3 for v = 0, 0.1, and 0.2, 
For the case of no central opening ( v s 0) the coefficient of Yj vanishes 
(Y|(0) i - “) and the solution becomes = 3j( Pji*") sin (t +0 ). Figure 4 
shows the pressure field (isobars) for the case v = 0 (no central opening), 
while Figure 5 shows isobars for v = 0.2. 

2 

Second Order Solutions . Setting the coefficient of e in Eq. (1) equal 

( 2 ) 

to zero yields the following equation for the second order potential <p 


( 12 ) 












t / 
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where only the known first order solutions appear on the right hand side. 
Following Maslen and Moore, Eq. (12) is assumed to have a solution of the 
form 

= A sin 2 (t + ne ) [r^ ( a) - (y + 1) f (a )] (13) 

Substituting Eq. (13) into the left-hand-sicic of Eq. (12) and the first order 
solution given by Eq. (7) into the right-hand>side of Eq. (12) yields thip 
following inhomogeneous Bessel equation for f ( a ): 


2 

f"(a)^ + 4 (l- S5)f(cv) = R2(o) 


where the primes indicate differentiation with respect to or . 

2 

Equating the coefficient of e in Eq. (2) to zero gives the following 
relation for in terms of the first and second order potentials and 




Substituting the first and second order potentials (i.e. Eqs. (7) and (13)) into 

( 2 ) 

Eq. (13) gives the solution for ' as: 


y[p2o(®)cos 2 (t + n6) + P 23 (Q^)J 


where 


Pjo(«) 


2 2 




ORIG!NAl PAGE 1$ 

OF POOR QUAL.ITY 

and P23(“> = ; [R>)]*-;R*(«)(i-p) < 1 « 

The second order contribution to the pressure is thus seen to consist of a 
double frequency component (i.e.| P 2 o^® ) 2(t + n6 )) and a steady 

component (i»e., Pjj( »)). 

In order to complete the second order solutions, the function f(cv) 

must be obtained by solving Eq. (14)* Since this is a second order differential 

equation, the solution for f(cy) will contain two arbitrary constants. These 

constants must be obtained by imiK>$ing the boundary conditions at the inner 

and outer radii; namely, that the pressure perturbations vanish there. This is 

not as straightforward as it at first might seem, because the relationship 

between the dimensionless radial coordinate a and the dimensional radial 
# 

coordinate r involves the angular frequency (u . Since the frequency is 
usually dependent upon amplitude for nonlinear systems, the application of 
the boundary conditions becomes complicated. Therefore the determination 
of the frequency dependence upon amplitude and the appiication of the 
boundary conditions will be deferred until after the third order solutions 
have been derived. The solution for f( » ) is given in Appendix A. 

3 

Third Order Solutions . Equating the coefficient of c in Eq. (1) to 

3 

zero yields the following equation governing the third order potential ; 


momM 

OF POOR QUALITY. 


+ (y ♦ (y -1) 

where only the known first and second order solutions appear on the right 
hand side. As in the Maslen and Moore analysis^ this equation is assumed to 
have a solution of the form: 

s F j(a ) sin 3 (t + n9 ) + ) sin (t + n6 ) (20) 

where the functions Fj(o) and F^(c^) satisfy the following inhomogeneous 
Bessel equations: 

= S'“> 

The inhomogeneous terms R|( ot) and ) appearing in Eqs. (21) and (22) 

are determined by substituting the first and second order potentials (i.e., 

Eqs. (7) and (13)) into Eq. (19), and they are given in Appendix B. 

(3) 

Also the third order pressure contribution P Was derived by 
applying the perturbation analysis to the nonlinear relation between pressure 
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♦ 


and velocity potential (Eq. (12)) to obtaini 

(23) 

Substituting the firsts second, and third order potentials into Eq. (23) gives 
(31 

the soiution for P as: 



y PjQ cos 3 (t + n 0 ) - yP 3|^ cos (t + n 0 ) 




where 


PjO^ O') = 3Fj + 




"s*'" > = >=5 - V 


(25) 


(26) 


The third order contribution to the pressure is thus seen to consist of a third 
harmonic component (i.e., P^q cos 3(t ^ n6)) and a correction to the 
fundamental (i.e., P^^ cos (t -f n 0)). 


■ 


In order to complete the third order solutions the radial functions 
F|(or ) and ) must be obtained by solving Eqs. (21) and ( 22 ) subject to 
the boundary conditions applied at the inner and outer radii. This depends on 
the amplitude dependence of the frequency which will be considered in the 
next section. The solutions for Fj( o) and F^( ot) are given in Appendix B. 

( 2 ) 

In both the second and third order solutions, the expressions for 9 , 

1 * 20 ^ P 3 ^(o') differential 

equations governing f(o), Fj(o) and F^^(o) agree with their counterparts 
in the original Maslen and Moore analysis if one replaces ) with 3^{oi ), 

Boundary Conditions and Dependence of Frequency Upon Amplitude. 
The perturbation technique of Maslen and Moore also requires that the 
frequency varies with amplitude according to the expansion 

Hi = 0)^ + e0)j+ u>2 + + (27) 

where u) ^ is the acoustic frequency. This is equivalent to expanding the 
eigenvalue | » bto /c^, thus 

p s Pq s Pj s P^ * P 3 ••• (28) 

In the remaining discussions the subscripts denoting the tangential and radial 
mode numbers n and k for the eigenvalue P will be dropped. 

In order to determine the coefficients P^, p 2 » •“ in Eq. (28), the 
pressure was required to satisfy the boundary condition P - 1 = 0 at the 


.. » > 
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outer radius o s p and at the inner radius o » v 9 f 

cP^^\p, 0 ,t)+ eV^Vp,e»t)-*- fV^^(p, 0 ,t)sO (29a) 

e P^^^(v P ,0 »t) 4. «M^^( vp ,0 ,t) + V p ,0 ,t) - 0 (29b) 

> 

Substituting the expansion of p given by Eq. (28) into Eqs. (29), expanding 
each term in a Taylor series about the acoustic value P^ or vp^, and 
equating the coefficients of each power of « to zero yield the boundary 
conditions for the second order function f(Q') and the third order functions 
F^(a) and 

In particular, equating the coefficient of e to zero yields R^( P ^) 

s R^( vp ^)=0 which is already satisfied from the definition of Py 

2 

From the coefficient of c one obtains the following expressions: 

Pj R},( p ^ cos (t + n 0 ) + PjqC Po^ 2(t f n 0 ) + P 23 (P ^ = 0 (30a) 

i- 

V Pj R^( V p^) cos (t + n 0 ) + P 2 q( vp^) cos 2 (t + n 0 ) + P 23 ( v P^) s 0 

(30b) 

Equations (30) must be satisfied for all values of t and 0 , thus each term 
must separately vanish* The first term is zero only if P^ vanishes, since 
tlva derivatives of the eigenfunctions R|^( P^) and Rp('^ Pq) are not zero at 
the boundaries (see Figure 3). The other two terms require that the second 

order pressure functions ^^23^°^^ vanish at = "^P^ and Qf s | 

\ 

P^. Thus the first order correction to the eigenvalue P| is zero which 
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agrees with the well-known result Irom linear acoustics that the frequency 
is indepen<lent of amplitudot From Eq* (17) the boundary conditions on 
P 2 q(o) yield the following boundary conditions to be satisfied by the second 
order function t(a )t 

o> " C ^ ]* 

i( nfj = 

On the other hand, Eq. (18) and the boundar/ conditions on require 

that R^( Pjj) and r'( vp^) vanish, which is impossible since the 
derivatives of the eigenfunctions R^ are finite at the boundaries. This 
implies that the condition that the pressure perturbation vanish at the inner 
and outer boundaries cannot be satisfied in the framework of the Maslen and 
Moore theory. 

The failure of the second order solutions to satlsiy the boundary 
conditions is the difficulty mentioned previously. This residual pressure 
perturbation (P - 1), which is independent of time and angle, is given to 
second order by (using Eqs. (29), (16) and (18))t 


P'O.) = 

(32a) 

p'(ve,)= -J*' [ 

(32b) 
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where P' Is the dimensionless pressure perturbation P-1* These boundary 
residuals appear to arise from the nonlinear terms in the relationship 
between pressure and velocity potential (i.e.^ Eq. (2))* It Is possible that 
these boundary residuals arc so inherent error associated with the Maslen 
and Moore analysis when one attempts to satisfy pressure boundary 
conditions with a perturbation analysis based on a velocity potentiai . In the 
original Maslen and Moore analysis tor a cylinder with rigid boundariesi one 
required the normal velocity component to vanish at the boundary* Since the 
velocity is linearly related to the velocity potentialr the boundary residuals 
did not appear in the original analysis. Further analysis of this problem has 
not yielded a method of satisfying the pressure boundary conditions to 
second orderi Therefore the approach to be taken in the remainder of the 
analysis will be to accept these boundary residuals as part of the error 
incurred by using an approximate perturbation analysis. 

Returning to the Taylor series expansion of Eqs. (29), setting the 
coefficients of e equal to zero yields: 

[ 2 **34^ ^ o^] + u6) + P 3 q( cos 3(t -t- n 6 ) s 0 

(33b) 

where ^ j = 0 has been used. Again each term must vanish separately. The 
second term, along with Eq. (25), yields the following boundary conditions to 
be satisfied by Fj( ah 


(34a) 


t 
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Fj( vp^y B rJ^(v P lU vp^) (JPb) 

The first term and Eq. (26) give the boundary conditions for P^(o) in terms 
of the second order correction to the eigenvalue P 2 asfoliowst 

Fj( P^) s Rj^( p^)f ' ( p^) - pjRj^f p^) (?5a) 

PjCvP^)** ^R'^(vp^)f'(vp^) - vP2R^(vp^) (35b) 

.The eigenvaiue correction p . can only be obtained after the 
inhomogeneous Bessel equation for F^(cr) (i«e.» Eq. (22)) has been solved. 
This analysis is given in Appendix B and yields the following expression for 
p2« 

“n<eo>- 

(36) 

where a a 3JPJ/JJ''PJ and F«(o) is the particular solution of Eq. 

non o 5p ^ ^ 

( 22 ) given in Appendix B. 

Total Solutions . In order to obtain numerical solutions relevant to the 
spinning wave engine, expressions for combining the first, second, and third 
order solutions given in the preceding sections are needed. These are based 
on the expansions of the velocity potential and pressure given by Eqs. (3) and 
(4) where the dependence of frequency upon amplitude is accounted for by 
the Taylor series expansion procedure used in the previous section. Since we 


23 


are interested in interior points as well as boundary points, the Tayioi’ 
expansion is performed about 3^r where » s(S^+ t 

The pressure perturbation p' is thus obtained from Eq, (4) using a 
Taylor expansion about p^r to givet 

P'(r,0,t)c eP^^\p^r,9,t) 4^ P^^\ p ^r, 6 ,t) 

+ P^r,e,t) ♦ Pg** ^ 

Using Eqs. (10), (16), and (24) in Eq. (37) yields 

P s - Y {e ^ + *^**34^ ^ o' C 2 0 ® ^ 

+ Pjo( <^) cos 2(t + n 6) ^30^“ ^ 3(t 4 n0 )} (38) 


where ^2 ~ ^2^^o computations are made using the acoustic 

eigenvalue P ^ in determining <y (i.e., using 01 = p^r in Eq. (38)). 

The acoustic velocity components are aiso of considerable interest in 
the spinning wave engine analysis, since they are needed to determine the 
acoustic pumping of a steady fiow through the engine. The radiai velocity v 
and tangential velocity w are determined from the definition of the 
velocity potential and Eq. (3) as follows^ 




boT 


( 2 ) 


a^(3) 

TqT 


(39) 


t 



2 

where again cv = ^ the Taylor series expansion 

about ^^r and neglecting terms higher than third order in < gives: 


v(a,e,t)= (Pj,r)+ V 5?”^ 




(3) 


«H) 




w(o,e,t) = { * 4|- < ^o'’^ “11“ ^ < V ^ 

.( 1 ) -. 2 ^ 1 ) 


r Sv' ' f B r> 1 Sr ^ 'O r)1‘l 


m 


Introducing expressions for the velocity potentials from Eqs. (7)t (13)» and 
(20), performing the indicated differentiations, and grouping like harmonics 
together yieids: 


vCa.e.t) = {« R'(a)+ [Fj(ff)- C,!! } ‘‘"(•♦'’S* 

♦ €*rR„(a)R'„(a ) - )j sin 2(t ne) 

+ F'j(a ) sin 3(t + ne ) (i>3) 


w(tf,e,t) * «R„(a) f Qji<3iR'^(a)~R^ia))]}coi(Une) 

^ 5 *^[ ^ ^ Y + Of ( Q^) J cos 2(t ■•■ n e ) 

F j( ot) cos 3(t n e ) 

These expressions can be used to determine the magnitude of the 
steady flow driven by the spinning wave. Since only the radial component 
contributes to a through flow, only Eq. (43) for V is of concern here. The 
steady through flow is obtained by taking the time average of v over one 
period of oscillation. Since the radial velocity consists entirely of terms 
which are sinusoidal in (t + n6) (fundamental with second and third 
harmonics), the time average of v is zero. There is no steady or time- 
independent term. Thus the Maslen and Moore analysis predicts no steady 
flow generated by nonlinearities in the wave motion when the pressure 
perturbation is required to vanish at the boundaries. This result appears to 
be independent of the difficulty in satisfying the boundary conditions for the 
second order solutions. 

4 

Typical Numerical Solutions . The expressions derived above for the 
pressure perturbation (Eq. (38)) and the velocity components v and w 

(Eqs, (43) and (44)) were used to obtain numerical solutions for a typical 
spinning wave engine configuration. In these calculations an annular pancake 
geometry was assumed with the ratio of central hole diameter to disk 
diameter v of 0.2. The working fluid was assumed to be air with a specific 
heat ratio v of 1.4, since combustion is not considered in this analysis. 
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Finally a first tangential (l.e., n = X> k » 1) spinning mode was assumed with 
an amplitude parameter c of 0k3i 

Before the acoustic pressure and velocity can be calculated the 
following functions are needed: the radial eigenfunction and its 

derivative, the second order function f(o ) and its derivative, and the third 
order functions Fj(a ) and F^(ot ) and their derivatives. The radial 
eigenfunction for the first tangential mode has already been shown in Fig. 3 
where 6^ = <>.2357. The second order functions f(a) and f^(cv) are shown 
in Fig. 6 for v = 0.2 and Y = 1.4. From these and the radial 
eigenfunctions, Eqs. (17) and (18) yield the second order pressure functions 
P2 o^Q') and P23(cv) which are plotted in Fig. 7. In order to obtain the third 
order functions Fj((v) and Fj(a), Eq. (36) was used to yield ^2 = 
-^49254. The function F^(a) and the corresponding pressure function 
P3o(cy) (obtained from Eq. (25)) are shown in Figure 8, while Figure 9 shows 
Fj(a) aod the function ) + a C2 (obtained from Eq. (26) and 

(38)). 

Pressure waveforms for the spinning first tangential mode were 
computed using Eq, (38) for two radial stations: r = 0,56 where R|(a ) has 

a maximum (pressure antinode) and r = 0.8. These waveforms are plotted in 

3 3 

Fig. 10 for both second order ( e terms neglected) and third order ( e 

terms retained) approximations with e = 0.3. At r = 0.56 the second order 

waveforms are only slightly distorted from a sinusoidal shape since P2 q(« ) is 

relatively small, while at r = 0.8 the second harmonic distortion is readily 

apparent for here P2 q(°^ ) attains nearly its maximum magnitude. Adding 
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Figure 9. Third Order Function F^Ca) and Corresponding Pressure 
Function fof Spinning IT Mode. 






the third order terms tends to f latten the peaks at r s 0.56 and smooth out 
the minima at r = 0.8. At the boundaries r = 0.2 and r s 1.0 the 
oscillatory component of the pressure vanishes, however a time independent 
component or mean pressure shift remains. This boundary residual which was 
discussed previously (Eqs. (32)) amounts to about - 2% of the ambient 
pressure at the inner boundary and about - 0.5% of the ambient pressure at 
the outer boundary for c = 0.3. 

The radial velocity waveforms were calculated using Eq. (43) at the 
boundaries where they attain their maximum amplitude. These waveforms 
are shown in Figure 11 for both second order and third order approximations. 
The principal effect of nonlinearities on the velocity waveforms is to 
steepen the rising portion of the curve and flatten the descending branch 
while preserving symmetry about zero velocity. This is more apparent at the 
outer boundary r = 1. Due to this symmetry, there is no steady flow through 
the system. 

The effect of amplitude upon frequency is determined from Eq. (27) 
as: 



2 ^2 
• i; 


or 



(45) 


For the present case, ^2 ~ “ *49254 for which C 2 = " *ii628. Thus 
frequency is seen to decrease as the amplitude parameter is increased. A 
similar variation of frequency with amplitude was calculated by Maslen and 
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Figure 11. Boundary Radial Velocity Waveforas for the Spinning IT Mode. 
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Moore for the hard>waUed cylinder with 

Irequency shift is quite smail^ giving a frequency only about one percent less 
than the acoustic frequency. 

Analysis with Steady Potentials 

The straightforward application of the perturbation technique of 
Masien and Moore to the open annular cavity does not yield an acoustically 
driven steady througti flow as shown by Eq. (43). In this analysis the solutions 
assumed for and were patterned after the solutions 

assumed by Masien and Moore which did not explicitly contain a steady flow 
term. In this problem* hov^ever* a steady flow term is permissible since the 
boundaries are open. In an attempt to predict acoustic pumping, the Masien 
and Moore approach will be modified by including a steady potential in the 
solutions for (p^^^and <p^^\ 

The analysis is modif ied by adding a steady potential to each of the 
unsteady potentials used previously. Thus the velocity potential becomes 

?<»]* 

(«) 

where the steady potentials are functions of o; and e and the sum of 
the potentials cp^^^ + must satisfy Eqs. (5), (12), and (19). The existence 
of V however, implies the existence of acoustic pumping in the first 
order or linear regime (small amplitude). Since acoustic pumping i$ expected 
to arise from the nonlinearities of the system, it will be postulated that the 
first order steady potential is zero. Thus the first order solution for the 
velocity potential is still given by Eq. (7) and the first order pressure is still 
given by Eq. (10). 


The second order solutions are determined by replacing in the 
ieft-hand-side of Eq. (12) with and using Eq. (7) in the right- 

hand-side for This yieids as given by Eqs. (13) and (14) while 
must satisfy Laplace's Equation) 


s 0 


(47) 


Due to symmetry considerations the steady fiow driven by a spinning 

acoustic mode is CKpected to be independent of the angular coordinate 6 . 

Soiutions of Laplace's Equation which satisfy this criterion are the potentiai 

source or sink which is a purely radial flow and the potential vortex which is 

a purely tangential flowt Although only the radial flow component can 

contribute to a mean flow through the annular cavity, there is no apparent 

physical reason to exclude the vortex flow component. Since the acoustic 

wave rotates around the chamber, it is quite likely that a vortex flow is 

generated. Since Laplace's Equation is linear^ superposition of solutions is 

.( 2 ) 

permitted, and the steady potential 9 ' will be assumed to be a linear 
combination of source and vortex flows given by; 

.-«>= (*») 

where K and K are constants to be determined. 

5 V 

( 2 ) 

The second order contribution to the pressire P ^ is obtained by 
substituting and into the right-handUside of Eq. (15). Since 

steady component and = 0 (by definition) the 

(2) _l2i 

expression for P is unaffected by the second order steady potential cp' • . 


* * 4 
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Thus P^^^ is still given by Eqs. (16), (17) and (18). From this result It is seen 

that including the second order steady potentials does not remedy the 

dif£iculty in satisfying the pressure boundary conditions to second order. 

The third order solutions are obtained by introducing Eq. (7) for 

and given by Eqs. (i3) and ((>8) into the right-hand*side of Eq. 

(i9). Expanding the right-hand-side, collecting terms, and using + 9 ^^^ 

(%\ 

in the left-hand-side yields the solution for 9 ast 


9 ^^^ s Fj(a ) sin 3(Une ) + F^(«) cos(t+ne ) + F^ia ) sin(t+n e) 

(^9) 

where Fj(a) satisfies Eq. (21) as before, but, F^( a) and F^( o) are solutions 
of the following inhomogeneous equations: 


f'.: + s F' + (l- !£)F.= ^ R'<«) 




a * 4 


a 


2 ' 4 ^ O' 


f '' ♦ ^ )p5= «5(«)- 


(30) 


(31) 


vfhere Rj(«) is given by Eq. (B-16) as before. The steady potential 9 ' ^ 
must also satisfy Laplace’s Equation: 

7 ^ 9 ^^^ »0 (32) 

The third order contribution to the pressure is obtained from Eq. (23) 
using Eqs. (7), (13), (98) and (99) to yield: 


s - Y I P3o^® ) c®* 0 ) F33(« ) sin(t-»-ne ) 


■f Pj^Cor) cos(t4tie); 


(33) 
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where ) is given by Eq. (25) as before^ while P 33 (°') Pj^( a) are 
given by: 

K. . 

Pj^((v) = -F^(a) + 0tf) 

nK 

P3^(a) = Gy(a) + (55) 

where G^(Qf) is the right-hand-side of Eq. (26). 

It is seen from the above analysis that including a steady term in the 
-(2) 

second order potential 9 ' affects only the third order contribution to the 
pressure through the new term «) sin(t+n 9 ) and by modifying 
This effect can be viewed as shifting the phase and changing the amplitude 
of the third order pressure component which oscillates at the fundamental 
frequency. A similar effect occurs in the third order potentiai through the 
new function F^( o) and a modification of F^(oi'). 

The boundary conditions are applied using the Taylor series 
expansions as before with the eigenvalue 3 expanded according to Eq. (28). 
Again the first order correction p j vanishes and the boundary conditions on 
f( «) are again specified by Eqs. (31). The third order function F|(o) also 
satisfies the same boundary conditions as derived previously^' namely, those 
given by Eqs. (3^), The boundary conditions on F^(q^) are also the same as 
before, since the contribution of the vortex flow to Pj^(o' ) vanishes at the 
boundaries (see Eq. (56)). Thus the boundary conditions for F^(o) are given 
by Eqs. (35) in terms of 62 * However the solution for F^( a) depends upon 
the vortex flow through the additional term on the right of Eq. (51). The 
pressure boundary conditions also require that ^33^ vPo^ **33^^ o^ "" ° 
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V I 


which Irom Eq, yields: 

V y = ^ “n' ^o> 

Solutions for F^(o) and F 5 (q') are given in Appendix C. It is also 
shown in Appendix C that the boundary conditions do not provide enough 
information to uniquely determine the frequettcy correction p 2 > steady 
flow parameters K and and the constants of integration appearing in 
the functions F^lo) and F^(a), Furthermore, it is shown that there is no 
value of K for which Eqs. (57) are satisfied. 

From the above analysis, It appears that adding steady flow terms to 
the Maslen and Moore analysis does not provide a means for predicting 
whether a large amplitude spinning wave can generate a steady flow through 
an open annular cavity. 



TEST MODEL ' ND INSTRUMENTATION 

Test Model 

An illustration of the spinning wave engine test model is in Figure 12i. 
The assembled model is shown in the photographs of Figure 13, Two parallel, 
61 cm (24") diameter discs are used to form the cylindrical cavity for the 
spinning wave. The cavity is open to the environment around the entire 
periphery and at central ports 12.2 cm (4.8") in diameter ( v = 0,2). The 
spacing between the discs is variable via the four threaded support rods 
located outside the cavity. One disc is fabricated from 2.5 cm thick, 
transparent plexigias to allow flow visualization. The other disc is fabricated 
from 1.25 cm thick aluminum plate. Four acoustic drivers, at 90° intervals, 
are mounted to the aluminum disc with the active surface flush with the 


internal surface. The model is equipped with ten ports for flush mounting 2,5 
cm microphones. Three of these, which are clearly visible in the 
photographs, are along a radial line on the plexigias disc. Seven are on the 
aluminum disc along a radial line and at three circumferential positions at 
90° intervals. The radial location of these circumferential microphone ports 
and the acoustic drivers coincides with the pressure antinode predicted by 
the acoustic solution. Also shown in the photographs, but excluded from the 
illustration of Figure 12 for the sake of clarity, are five static pressure taps 
along a radial line on the plexigias disc. 

Figure 14 presents a plan view of the apparatus and shows locations 
of the various instrumentation and the corresponding nomenclature used 


herein. The four drivers are indicated as 


Dl^, Dig, D2,, and D2g. The 




PPI 



Figure 12. Spinning Wave Engine Test Model. 
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Figure 1 3. Test Model Assembly. 






microphones, Ml, M2^ and M3, are shown at the primary locations for the 
spinning wave tests. The static pressure ports indicated hy Pj through are 
equally spaced along a radial line. P^ is located on the surface opposite to 
that of D2|^. Finally, six locations >used for the constant temperature hot 
wire measurements are shown as W1 through W6. These cover one quadrant 
of the flow field with spacings at intervals at both the Inner and outer 
radius of the cavity, 

Instrumentation 

Four 75 watt University drivers. Model ID-75, were used to excite the 
oscillatory flows. In most cases, the four drivers were operated in pairs with 
Dly^ and Dig as one pair and D2^ and D2g as the other. Each pair was 
powered by a single Krohn-Hite 50 watt amplifier. Model DCA-50R, and 50 
watt matching transformer, Model MT-56. The two drivers of a pair were 
excitj^d 180° out of phase by reversing the polarity of input leads. A 
Hewlett Packard Model 203A two-channel, variable phase function generator 
was used to provide simple harmonic input signals to the two amplifiers. 
Operation with any combination of from one to four dHvers was possible by 
connecting or unconnecting input leads. 

Three B & K condenser microphone cartridge type 4134 were used for 
measuring the oscillatory pressures. The microphones could be relocated to 
any of the ten microphone ports shown in Figures 12 and 13 by interchanging 
microphones and cylindrical plugs. Microphone outputs were measured by a B 
& K type 2606 amplifier with frequency range from 2 to 200,000 Hz. 
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Hot-wire measurements were made wit^ a TSl single channel hot- 
wire anemometer system consisting ot a constant temperature anemometer 
module, a monitor and power supply module, and a signal linearizer with a 
4th degree polynomial curve fit. A tungsten platinum coated hot wire with a 
sensor diameter of , sensor length of 1.25 mm, and upper frequency 
response of better than 200 kHz was used. Hot-wire outputs were measured 
both with an averaging voltmeter and a true RMS voltmeter. 

Input .and output AC signal wave forms were observed and 
photographed using a Tektronix type 549 oscilloscope equipped with a 
Polaroid camera attachment. Furthermore, phases between the various AC 
signals were measured with a Wavetek phase meter Model 740. 

Mean or DC static pressures were measured using a highly stable and 
accurate Datametrics Barocel Electronic Manometer with a 10 mm Hg 
variable capacity pressure transducer. Measurements are accurate to 1% of 
reading over the range from about 95 to 175 dB. 

Calibrations 

The microphones were calibrated with the standard B & K 
Pistonphone type 4220. The hot-wire calibration of voltage output versus 
flow velocity was provided by TSl. This calibration, which was made using 
their low velocity, steady flow calibrator, covered the range of from 0 to 6 
m/s. The calibration data was fitted to a fourth degree polynomial of 
velocity versus voltage within a velocity accuracy of better thah 196 for all 
test conditions reported herein. The analog circuitry of the anemometer 
linearizer were set using the coefficients of the polynomial so that the 


output voltage was equal to the velocity. Thus, assuming quasi>steady flow 
over the wire, as is the practice for turhulent flow measurements, the 
linearlzer output yields a faithful represehtatiofi of the velocities and 
velocity wave forms. 

Since the drivers are mechanical devices with inertia the output and 
input signals cannot be perfectly in phase. Of more importance, differences 
in the mechanical drivers will cause differences between the phases of the 
output and input signals. Checks of the drivers used in this investigation 
indicated that the output signals from the various drivers could be out of 
phase by as much as 5^ for the same input signal. This is relatively small in 
comparison with the 90^ and 180^ phase shifts used for the tests reported 
herein. Nevertheless, it does mean that it was not possible to obtain a 
perfect spinning wave, but, instead, the superpositioning of a predominant 
spinning wave on a weak standing wave was obtained. However, the 
existence of a superimposed weak standing wave does not alter the meaning 
or conclusions from the tests. 
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EXPERIMENTAL RESULTS AND DISCUSSION 
The experimental investigations include measurements of the 
resonant frequencies within the cylindrical cavity, phase changes for 
transition from standing to spinning wave modes, radial distributions of the 
osciiiatory (AC) and mean (DC) pressures, and simultaneous hot-wire and 
microphone data. In addition, flow visualization studies were made as a 
qualitative verification of the hot-wire results. The results are presented 
and discussed in the following sections. It is important to point out that 
exploratory studies yielding meaningful qualitative rather than more exact 
quantitative results were emphasized in order to maximize the benefits for 
the time available. 

Resonant Frequencies and Transition from Standing to Spinning Waves 

Natural frequencies of oscillations in the cylindrical cavity are 
indicated by the results of Figure 15 which presents plots of the RMS 
pressures measured by two microphones as a function of driver frequency. 
Results are shown for two power levels. For each power !evel» the driver 
input power was maintained while the excitation frequency was varied. 
Resonance occirs where the RMS pressures peak. 

Figure 15a Shows results for the standing-wave mode of operation 
excited by the single driver Dl^ (see Figure 14 for driver and microphone 
nomenclature and locations). The first harmonic occurs at a frequency of 
692 Hz. The second and third harmonics are at frequencies of 1065 Hz and 
1275 Hz, respectively. 
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Figure l^b shows results for the spinning wave mode of operation. 
Data are included for the spinning wave excited both by the two drivers 01^ 
and D2^ and by all lour dHvers (Dl^i Di^t and D2g). In both cases> 
02^ is set to lag Dl^ by 90^. Also, when using the four drivers, D1 g lags 
Dl^ by 180® and D2g lags D2^ by 180® since, as discussed before, the input 
leads for these pairs of drivers are reversed. The resonant frequencies are 
the same a$ those of Figure 15a for the standing waves. This Is to be 
expected since these spinning waves represent the superpositioning of the 
two standing waves oscillating 90® out of phase. 

Additional results verifying the natural frequencies are presented in 
Figure 16. This flgire shows power spectra computed from the outputs of 
microphones Ml and M2 when white noise is used to excite driver D2^. 
These results were obtained by feeding the real time microphone outputs 
directly to a Hewlett Packard Fourier analyzer and are based on an 
ensemble of 100 data strips and a frequency band width of 20 Hz. Because of 
a slight error in the plotting scale factor the frequencies are shifted slightly 
toward ihe higher values (e.g., the cutoff frequency should be precisely 
2500 Hz rather than the indicated value of about 2575 Hz). Accounting for 
this slight shift, the resonance frequencies for the first three harmonics, as 
indicated by the peaks in the power spectra, are about 700, 1070, and 1275 
Hz, in excellent agreement with the results of Figire 15. Also plotted on 
Figrre 16 are computed curves of the magnitude and phase of the transfer 
function which relates the outputs of Ml and M2. These results show that 
the two microphones respond in phase and with nearly a one-to-one 
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relationship (the magnitude of the transfer function is slightly 9 *eater than 
1.0 due to differences in the microphones). This and the similarities in the 
pov/er spectra attest to the quality and symmetries of the model and 
instrumentation. 

The experimental value of 692 Hz for the fundamental frequency is 
significantly different from the theoretical acoustic value of 760 Hz for v s 
0.2. This cannot be attributed to higher order effects. At the highest 
pressure level for resonance (RMS pressure s 157 dB) the value of the 
amplitude parameter e , as used in the perturbation analysis, is only 0.025. 
At the lowest pressure level for resonance at the fundamental frequency the 
value of the amplitude parameter is less than 0.01. Thus, second and higher 
order terms must be small. It is noted, however, that the experimental 
frequency is nearly the same as the theoretical acoustic value for a model of 
the same outer diameter but without a central hole (i.e., v s 0). This might 
indicate that the constant-pressure boundary condition is not satisfied at the 
inner boundary. 

Figure 17 is included to illustrate phase changes within the cavity for 
transition from standing to spinning wave modes of operation. These results 
were obtained using drivers Dl^ and only, and with D2^ lagging Dl^ 
by 90^. Similar results were obtained using all four drivers but the results 
were not completely documented. The phase lag between the output of Ml 
and M2 is plotted against the ratio of the input power to D2y^ to that of Dl^. 
The power input to 01^ was held constant while the power input to D2^ was 
varied. With a power ratio of zero a standing wave exists and the theoretical 





phase lag is 180*^. As the power ratio increases two standing waves 90® ap&rt 
and 90® out oJ! phase are superimposed to produce combined spinning aind 
standing waves. For a power ratio of unity a pure spinning wave is obtained 
and the theoretical phase between the two microphones, is 90^, 
corresponding to the angular displacement of 90®. Also shown on Figure 17 is 
the theoretical curve for the phase lag derived by superimposing standing 
waves. It can be shown that this theoretical phase lag is the arc-cotangent 
of the power ratio. The agreement between the theoretical and experimental 
results shown on Figure 17 is good. 

Radial Distributions of Pressures 

Figure 18 presents radial distributions of the RMS value of the 
periodic component of pressure for the spinning wave mode of operation. 
The pressures are normalized by the value at r = 1.0 in order to facilitate a 
comparison with the theoretical solution. Included in Figure 18 is the 
acoustic solution which is appropriate for these experimental data since the 
maximum value of the amplitude parameter e is less than O.ul. As shown 
by Figure 18, the experimental pressures are in good agreement with the 
acoustic solution for 0,33 < r< 0.87 and it seems that the periodic pressure 
will decay to zero at the open boundaries where r = 0.2 and 1.0. However, 
the constant pressure boundary condition is not adequately confirmed, and, 
in retrospect, the test model should have been equipped with microphone 
ports nearer these open boundaries. 
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Figure 19 presents radial distributions of the diffek*ence between the 
mean or steady component of pressure, p , and the ambient press*jre, p^. 
These time-averaged press^ire differentials were measured using a slow- 
response (i«e., high volume) pressure transducer. The five data points in the 
central region of the cavity from r ? 0,35 to r = 0.85 were measured using 
the static pressure taps labeled P. through in Figure 14. The two data 
points indicated by the triangular symbols were measured at the open 
boundaries ( r s 0.2 and r = 1.0) using a 0.15 cm static pressure probe. 
These are the only data that indicate negative pressures. The accuracy of 
these two data points is questionable since crossflow over the probe due 
either to oscillatory motion or boundary layer streaming would reduce the 
pressure readings. Thus, the actual pressure at these open boundaries may be 
slightly higher than that indicated by these data points. Regardless, the 
mean pressures at the boundaries are nearly ambient. 

In the upper portion of Figire 19 the pressure differentials are 
normalized by the amplitude of the oscillatory pressure, p' , measured at r 
= 0.6 (i.e., nearly the maximum value of the oscillatory pressure). In the 
lower portion the pressure differentials are normalized by the ambient 
pressure. Results are shown for the spinning wave <friven both by 4 and 3 
drivers. The reason for including the case of 3 drivers is that driver D2^ is 
located along the radial line of the static pressure taps on the disc opposite 
these taps. It was thought that this driver might have a direct effect upon 
the pressures and, thus, driver D2^ was not used in the tests with 3 drivers. 
As shown by the results of Figure 19, driver D2g does significantly change 
the pressure magnitudes, but it does not change the qualitative nature of the 
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data. 

There are two important points to be made concerning the results of 
Figure 19. Firsts the oscillatory motion associated with the spinning wave is 
accompanied by a radial distribution of mean pressure. The theory also 
predicts a steady component of pressure superimposed on the oscillatory 
pressure. The magnitudes of the theoretical values of the steady component 
of pressure (i.e„ (p - P^VPq = - Y P 23 ®^) fhe boundaries at r = 0.2 and r 
s 1.0 are about the same as the maximum value within the cavity (see Figure 
7). Also the theory predicts that the steady component of pressure is greater 
than ambient near the mean radius, less than ambient at the open 
boundaries, and equal to ambient at about r = 0.38 and r = 0.78. However, 
in contrast, the experimental results show that the mean pressure is greater 
than ambient for the entire cavity except perhaps very near the open 
boundaries where it may be slightly less than ambient. The second point is 
that the pressure differential P - Pq is small in comparison with the 
amplitude of the oscillatory pressure p' , As shown by the upper portion of 
Figure 19, the ratio of these pressures is of the order of 10 . The value of 
the theoretical amplitude parameter t for these tests was of the order of 
10"^ also. Therefore, since p' arises as a first order term, P- P^ must 
enter as a second order term. This is as predicted by the theory which shows 
that the steady component of pressure first arises as a second order term. 
The amplitude parameter based on p^ = 161.3 dB for the tests with 3 drivers 
is about e = 0.03. Using this in the second order theory (i.e., in the equation 
(p " ^ ^23 ® ' yiel<ls a maximum theoretical value of (p ^ 

Po)/Po = 0.9 X 10"^. This is in fair a^-eement, but slightly lower than, the 
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maximum value of 1.25 x 10~^ for the corresponding results in the lower 
portion of Figire 19. 

Hot Wire Results 

The response of a hot wire in oscillatory flows and the experimental 

methods used In obtaining the hot-wire data will be considered first to aid in 

understanding and interpreting the hot-wire results. Figure 20 shows 

idealized responses of a hot wire in oscillatory flows. Figure 20a presents 

the case for a pure sinusoidal motion with zero mean flow velocity. The flow 

velocity u for one cycle of motion is shown in the upper plot. The lower 

plot shows the corresponding velocity u^ that would be measured by a hot 

wire assuming a faithful, quasi-steady response to the flow. For all results 

reported herein, this represents the idealized output voltage since the 

anemometer signals were linearized and scaled by the analogue circuitry so 

that the output signal was equal to the velocity in accordance with a steady 

flow calibration. Because the hot wire is insensitive to flow direction the hot 

wire yields a rectification of the actual velocity cirve, as shown in the 

figure. This gives rise to an apparent mean velocity u which can be 

w 

evaluated as the average or DC component of the output voltage. 

Figure 20b presents the case in which a mean flow velocity u is 

superimposed upon the oscillatory velocity. The rectified signal from the hot 

wire is altered so that the maximum velocity for the two half-waves, u^j 

and u ^2 > half-wave lengths, \ ^ and \ different. 

Furthermore, the indicated mean velocity u and the RMS value of the AC 

w 

component of the output signal change with u . For this case the mean flow 





velocity is given by 


u = (U„, - w„2>/2 

Also, it can be shown that u is determined by and the RMS value of the 
AC compnent of the output signal. As will be shown later, the hot wire does 
not reproduce the rectified velocity curve. Apparently the quasi-steady 
condition is not sufficiently Valid and, thus, the output signal is distorted, 
particularly during the time priod near flow reversal when flow conditions 
change most rapidly. Nevertheless, asymmetries in the output signal like 
that shown in FigU'e 20b is a useful measure of the existence of a mean 
flow. 

The asymmetries shown in Figure 20b reveal the existence of a mean 
flow velocity but not the mean flow direction. For the tests reprted herein 
the mean flow direction is determined by the phases between the hot wire 
and microphone signals. This is illustrated in Figure 21 which considers the 
signals from microphone M2 and the hot wire at location Wl. For this 
presentation it is assumed that the wave spins clockwise in the orientation 
of Figure Also note that Wl is located 90® clockwise from M2. Figure 
21a shows the time-dependent output voltage E for M2. This represents a 
reference trace for the lower plots. Figure 21b shows the pressures at M2 
corresponding to the microphone output trace. Because of the electronics 
the pressiM’e is 180® out of phase with the microphone output. The radial 
compnent of velocity, u_, at the inner opn boundary ( r = a) and at the 
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a) Voltage output from microphone M2 



b) Alternating pressure component at M2 


c) Radial velocity at r and at angular position of M2 
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d) Radial velocity at r and at angular position of W1 
^ u /~ “"2 



e) Idealized response of hot wire at location W1 


Figure 21. Illustration of the Relation between Microphone 
Output and Hot-Wire Response. 
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angular position of M2 is shown in Figire 21c. This velocity wave as shown 
leads the pressure wave of Figure 21b in accordance with theory. For the 
sake of iilustrationi it is assumed also that the oscillatory velocity is 
superimposed upon a positive mean flow velocity (i.e., there is a net inflow 
rate at r = a). Since the wire location W1 is displaced 90° in the spinning 
direction from Ml, the oscillatory velocity u^ at W1 lags that of Figure 
21c by 90°. Figure 21d shows u at location Wl. Finaliy, Figire 21e 
presents the rectification of this velocity at Wi and, thus, the idealized 
response of t'lie hot wire* The important point of this iilustration is the 
relation between the measurable hot wire and microphone signals which can 
be displayed simultaneously on an oscilloscope. As shown for this case of 
inflow at r = a the phase lag between the maximum veiocity u^^ and the 
maximum microphone output is 139° (3n/2). The corresponding phase 
lags for the hot wire located at W2 or W3 are 180° and 229°, respectively. 
On the other hand, if the mean flow velocity u_ was negative each of these 
phase angles would change by 180° since the locations of u^| and u^^2 
would be interchanged. Thus, the phase between and u^^^ can be used 
to establish whether u^ is positive or negative. 

A typical set of oscilloscope traces from the spinning wave tests is 
shown in Figure 22. Operation is with all four drivers at the fundamental 
frequency of 693 Hz. Because the camera was set for a shutter speed of 0.1 
sec each trace shown on these photographs represent the superpositioning of 
about 70 sweeps across the oscilloscope trace. Since the band widths of the 
superimposed traces are small the signals are stationary* Figure 22a provides 












the zero voltage references for the two channels. Figure 22b shows the 

simultaneous traces of the output signals from microphones M2 and M3 

which have zero reference traces 02 and 03| respectively. This shows that 

the wave spins clockwise from M2 to M3 with the proper phase lag of 90^ 

from M2 to M3. Figure 22c shows the simultaneous traces for M2 and the 

hot wire located at Wl. This corresponds to the case depicted in Figure 21. 

The zero references for M2 and Wl are 02 and 03| respectively. Note first 

that the hot wire signals are always positive and, thus, present a 

rectification of the flow velocities. However, the minimum indicated 

velocity is not zero and, in contrast with the idealized response shown in 

Figure 21e; the valleys in the traces are well rounded. This is attributed to 

the failure of the quasi-steady flow calibration for these regions where the 

velocities change most rapidly and flow reversal ocars (during flow 

reversal the wake flows back across the wire). The important point from 

these results is that the successive peaks in the hot-wire traces are 

different, indicating a mean flow velocity superimposed on the oscillatory 

flow. Furthermore, the maximum voltage signal from the hot wire lags the 

maximum value of the microphone signal by about 133° indicating a net 

inflow as discussed in the previous paragraph. For reference, the locations 

for u. and u_ based on the microphone trace are shown on the 
®max °mln 
figire. 

Figure 23 presents simultaneous traces of the output signals from M2 
and the hot wire for all six hot wire locations. The zero voltage references 
are the same as in Figure 22a. The upper row of traces is for the hot wire on 



Microphone ^id Hot-Wire Traces at 6 Hot-Wire 
Locations, f = 693 Hz. 
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the inner boundary r = a and the lower row presents corresponding traces 

for the outer boundary r s b. Figure 23a is a repeat of Flgire 22c, The 

locations of u. ana u^ and of Uu and Ul are displaced in 
“max “min max roin 

increments of 45° corresponding to the 45° increments between hot-wire 

locations, as discussed earlier. As shown by the upper row of traces the hot- 
wire output is in phase with the radial velocity component and there is a net 
inflow rate at each location on the inner open boundary. The riidial velocity 
components at the outer open boundary is 180° out of phase with that at the 
inner boundary and, thus, the locations of the maximum and minimum 
velocities for the lower row of traces are interchanged with respect to those 


of the upper rows. But the hot-wire signals have also changed so that the 

f 

maximum and minimum hot-wire outputs occur around u, and u^ , i 

°max °min j 

respectively. Since Uc corresponds to outflow these traces show that 

°max I 

there is a net mass outflow rate at the outer open boundary. The importam: I 


point is that these results are consistent and show that flow is pumped 
through the cavity by the spinning wave and that the mean flow is radially 
outward from the inner open boundary to the outer open boundary. 

Additional hot-wire measurements were obtained with the standing 
wave mode of operation dHven using DIA and D2B (180° out of phase with 
DIA). These results were simply recorded from the oscilloscope screen and 
not photographically documented. The results were as expected and 
consistent with those for the spinning wave mode. Along the radial line 
through DIA the hot- Wire signals indicated a mean flow from the inner 
boundary at W1 to the outer boundary at W4, The mean flow rate 



decreased as the wire location was rotated around the periphery to W2 and 
W5. Finally, along the radial line normal to that through the drivers (i.e., at 
W3 and W6), where there must be a pressure nodCf there was no measurable 
hot-wire signal. 

R,esults extracted from the photographs of Figure 23 were used to 
compute the maximum and minimum flow velocities, u^| and u^ 2 > 
the mean flow velocity as given by 


u = (u^, j - 


The results are based on the steady flow calibration of velocity versus 
voltage output and, thus, assume quasi-steady flow over the wire during the 
measurements of u^ and u^^ . The velocities were obtained by measiring 
the displacements of the hot-wire traces from the zero reference traces and 
then applying a calibration factor for the oscilloscope deflection which was 
acquired from a separate test. The results are tabulated in part a) of Table 
I. For these results the wire was located midway between the upper and 
lower disc at y = 0 (see the sketch at the bottom of Table 1) and the mean 
of the RMS pressures measured by M2 and M3 was 155.3 dB. This was 
about the highest possible power level without driver induced distortions. As 
shown in the table, the average values of u at the inner and outer open 
boundaries are 0.34 and 0.13 m/s. These are very small for this low power 
level. Each of these is about 15% of the corresponding amplitude of the 
oscillatory velocity evaluated as (u^^ + U^^ )/2 and also listed in the table. 
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Table 1. Hot-Wire Results, £ « 693 Hz, 
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The amplitude parameter e is about 0.02 for these tests at the RMS 
pressure of 155,3 dB. This indicates that the acoustic solution should be a 
good approximation for these test conditions. However, the amplitudes of 
the oscillatory velocity at the inner and outer open boundaries as predicted 
by the acoustic solution for 155,3 dB at the pressure antinode are 6.0 and 3.1 
m/s. These are about a factor of three higher than the mean values of 2.2 
and 0.91 m/s listed in the table. Also, it is important to point out that the 
ratio of u at the inner and outer boundaries is about 2.5 whereas the flow 
area ratio is 5. These differences cannot be explained, but perhaps are due 
to a combination of things such as inadequacy of the steady flow hot wire 
calibration, asymmetries, and viscous effects. Nevertheless, it is felt that 
the qualitative behavior is correctly represented. The mean velocity and 
the RMS value of u^ are listed in Table la also since they were evaluated 
separately from the oscilloscope data using an integrating voltmeter and a 
true RMS meter. These data are entirely consistent with those obtained 
from the oscilloscope traces. 

Part b) of Table 2 tabulates results for the hot wire at W 1 and four 
values of y. The corresponding photographs of the oscilloscope traces are 
shown in Figure 23c, for y = 0, and Figure 2<F, for y = + 0.5 and y = 0.9 cm. 
For the three locations within the 1.25 cm gap (i.e,, y = 0 and y = + 0,5 cm) 
the velocities are ssentialiy the same. At y = + 0.5 cm the band width of 
the superimposed Jaces has inoreased very slightly relative to that at y = 0. 
At y - 0,9 the wire is located outside the gap and the oscillatory signal is not 
stationary as evidenced by the multiple traces. However, it appears that on 
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the average the peaks and vaileys in the traces are the same and, thus, there 
is oscillation with no mean flow. 

Flow Visualization 

Flow visualization was employed to support the hot wire predictions 
of a mean flow radially outward through the cavity. Srnoke and a tufted 
wand were used. A very sensitive tufted wand constructed with very fine 
strands of fiberglas showed without question that ther<$ was inflow around 
the complete periphery of the inner open boundary. Outflow was indicated at 
the outer boundary but, because of the very low velocities superimposed on 
the oscillatory motion, these results were not so conclusive. 

It was not possible to use a steady flow smoke trail to visualize the 
flow since the velocity of the smoke source was significant compared with 
the mean flow velocity. However, prior to starting oscillations it was 
possible to fill the cavity with smoke which becomes stationary and diffuses 
very slowly. Thus, the cavity could be filled with smoke and, after a short 
settling time, the driver quickly excited to produce the spinning wave. Tests 
of this type showed that the smoke was rapidly exhausted through the outer 
open boundary once the spinning wave was excited. Thus again the existence 
of a mean flow from the center outward was verified. 


CONCLUDING REMARKS 


In order to assess the feasibility of the spinning Wave engine concept, 
the acoustics of an annular cavity was investigated using a theoretical 
analysis and an experimental model. The primary objective of the 
investigation was to determine whether a spinning wave excited in the 
cavity could pump flow through the cavity, which was bounded on the sides 
by two parallel discs and open to the environment at the inner and outer 
peripheries. 

The theoretical analysis developed under this project was an 
extension of the nonlinear analysis of Maslen and Moore to study the case of 
an annular pancake cavity with open boundaries. This analysis Involves a 
perturbation scheme In which the velocity potential and pressure are 
expanded in terms of an amplitude parameter. Solutions were obtained to 
third order in this amplitude parameter. These solutions yielded a steady 
contribution to the pressure perturbation which was second order in the 
amplitude parameter and did not vanish at the inner and outer boundaries as 
required by the open boundary conditions. Fur thermore the analysis did not 
predict a corresponding steady velocity component; that is, the analysis 
could not predict the acoustic pumping of flow through the cavity. 

The spinning wave engine concept has been tested experimentally 
using acoustic drivers to excite oscillatory motion in an annular cavity. The 
spinning wave was achieved by superimposing orthogonal standing waves 
driven 90° out of phase. The major thrust of the experiments was to 
determine if the spinning wave pumped flow through the cavity. Also some 
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tests were made to assess the theory. 

Simultaneous measurements of hot wire and microphone responses 
have shown that the spinning wave pumps a mean flow radially outward 
through the cavity. This conclusion has been verified by flow visualization 
using both a sensitive tufted wand and smoke as flow tracers. Measured 
radial distributions of the amplitude of periodic pressures were in good 
agreement with the theoretical predictions. Also, the experiments showed 
that the spinning wave gives rise to a radial distribution of mean pressure 
which is of second order, as predicted by the theory. In contrast with the 
theory, however, the measurements indicated that these mean pressures 
were positive throughout the entire cavity. Furthermore, the experimental 
value of the fundamental frequency was about 10 percent below that 
predicted by the theory. Perhaps these discrepancies between theory and 
test results are associated with the difficulties encountered in satisfying the 
boundary conditions with the theory. 

To further evaluate the potential of the spinning wave engine, tests 
at higher power levels are needed. Because of driver limitations, the results 
reported herein were essentially at the acoustic level. In fact, the maximum 
value of the amplitude parameter was about 0.025. Furthermore, the 
accuracy of mean flow velocities measured by hot wires in such a flow field 
which is dominated by oscillatory motion is open to question. Thus, for 
accurate mean velocity data additional hot wire calibration tests are needed 
or the velocities must be measured by a more reliable instrument such as a 
laser velocimeter. 


In spite of the limitations of the theoretical analysis and the 
experimental measurements, the primary objective of this research project 
has been met. The experimental detection of a radial mean flow driven by 
the spinning wave at relatively low acoustic amplitudes has demonstrated 
the feasibility of the spinning wave engine concept. It is anticipated that 
this acoustic pumping will increase as the square of the acoustic amplitude 
giving much larger mean flow velocities at the larger acoustic amplitudes 
expected in a spinning wave engine. Thus further development and testing of 
the spinning wave engine concept is definitely warranted. 
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Appendix A 

Solution for Second Order Function f(a ) 


The second order function f(a) is described by the inhomogeneous 
second order differential equation: 


l"<a)« il'(o) + )«<»)= Rhct) 

a ^2 n 


(A-1) 


where ) is given by Eq. (7). Inspection of this equation reveals that the 
change of variable z s 2ar will transform this equation to a Bessel equation. 
Thus Eq. (A-1) becomes: 





f = rJ(z/2) 


(A-2) 


The homogeneous version of Eq. (A-2) has the general solution: 


= C^3^^z) . C2Yj„b) 


which Implies that the homogeneous solution of the original equation (Eq. 
(A-1)) is given by: 

f,^(a)= Ci32n(2a) + C2Y2n(2a') (A-3) 

The complete solution is the sum of the solution of the homogeneous 
equation f j^( o') given by Eq. (A- 3) and the particular solution fp ( a). Thus 
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Thus the general solution of Eq. (A-1) can be written as followsi^^^ 


Ha) - 


[=i -/ 


Y2n<2=^> 
W<a) 


] 




<=2*f 


J2„(2^^) Rn<«) 

W(o) 


] 


(A-«) 


where W(a ) is the Wronskian given by 

W(a)= 2J2„(2cy ) Yj^<2a ) - 2y2„(2a )a^„(2a ) (A-5) 

and primes denote differentiation with respect to the argument. 


The indefinite integrals in Eq. (A-^) may be replaced by definite 
integrals by making a suitable adjustment In the constants and C 2 • 
Choosing the lower limit of integration at the inner boundary as v 3^, the 
general solution becomes: 


i{a) = 


where 

K|(o) 


K^ia) 


[C| - K2< 0 ^) j l2rS2a ) + ['C 2 + Kj(a)j Y2„(2or ) 




(A-6) 

a 

f 

J 2 n^ 2 x) R^(x) dx 

(A-7) 


W(x) 

7 

Rn(x) dx 

(A-8) 

'*0 

W(x) 
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The unknown constant' and Cj are determined from the 
boundary conditions imposed upon f(o) as given by Eqs. (31). Appiyingthe 
boundary condition at or s givesr 

'*■ ^2^2n^^^^o^ " 4( y+0 ^ 


since Kj( vP^) s K 2 ( « 0 from Eqs, (A-7) and (A-8), The boundary 

condition at « = gives: 

♦ K,<P,Wj„(2e„)-K,(e„)Y^(2P,) 

(A-10) 


Soiving Eq. (A-9) for C 2 gives: 


4fer)K<^Po>l - 


^2n<2%> 


(A-ll) 


Substituting Eq. (A-ll) into Eq. (A- 10) and solving for Cj then yields: 




']/[' 


- p2n<2Po> ‘ ‘'2n<2v^,)V2„(2^^) 


'] 


(A-12) 


#>ll 
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yp,) . K,(p.)Y^ee„).Kj(s„)3,„(2e,) 

For given volues of Y » v and n (tangential mode number), the eigenvalue 
is readily determined and the quantities appearing in Eqs. (A-11) and 
(A- 12) can be calculated. Equation (A- 12) is first evaluated to obtain 
which is then substituted into Eq* (A-11) to obtain C 2 > For v s 0.2, Y s 1*4 
and s 4.23575, the constants of integration become Cj » - 0.36495 and 
Cj- - 0.04275. 

Once the constants of integration are obtained, Eq. (A-6) is used to 
calculate £(» ). In this process, & numerical integration technique must be 
used to evaluate K|(o) and K 2 (" )« 

Irv computing the third order quantities, the derivative f^( o) is also 
needed. Differentiating Eq. (A-6) gives: 

f V) « 2 [C^ - K2(a) 1 J2n^2») + 2 [G 2 + K^( a)l Tj^a) 

-K^(Q')J2n(2a) + k J(«)Y2„(20 ) (A-13) 

where kJ(q') and K 2 ( o) are simply the integrands in Eqs. (A-7) and (A-8). 
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Appendix B 

Solution for Third Order Functions Fj and F^ 

Solution for F|(a) 

The third order function Fj(a) Is described by the Inhomogeneous 
second order differential equation: 


vhere 

* S “n<>' 

■ ^ 12( y+1) ^ + 6( Y -l)j R„f - 6 (y +l)R'„f 'I 


(B-1) 


(B-2) 


The change of variable z s 3o will transform Eq. (B>1) into an 
inhomogeneous Bessel equation given by 


£F, ,dF 

dz2 


The homogeneous version of Eq. (B-3) has the general solution given by: 


F^h<z) = + CjjYj^fz) 
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and returning to the original variable a yields the homogeneous solution of 
the original equation as; 


''lh<''' = <=„J3„(3a).C,3Y3„(3«) (W) 

Adding to this the particular solution yields the following solution 

of the inhomogeneous equation^^^ (Eq. (B-2)); 


Fi(a) J3n(3») 



f Y3j^(3o)R^(a)do “I 

L " J 

f w]TF) ” J 


+ Y3n'3o') |^Cj2 * J Wj^ST) "J 


(B-5) 


where WiCor ) is the Wronskian given by; 

Wj( o ) = 333^(3 o )Y3^(3 o ) - 3Y3^(3ff )J^^(3cy ) (B-6) 

and primes denote differentiation with respect to the argument. 


Replacing the indefinite integrals in Eq. (B*5) with definite integrals 
and adjusting the constants of Integration gives; 


F,(a> = [c,i-Iu(»)] 33n(3») 


• i. (Mh. T ( 
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where 


i„(“) 




a 

f Rj(x)dx 

J W,(x) 

(B-8) 

a 

r Yj^(3x>R^(x)dx 

J -~w-:t5) — 

(B-9) 



The constants Cjj and C |2 are obtained by applying the boundary 
conditions imposed upon Fj(o) aS given by Eqs. (3<»). This yields the 
following equations to be solved for j and Cjj; 

‘^1 l’3n<’'' ^o> * ^12’^3n<^ '' = - i? < V *l)Rn('> P Jl' ( '■ (B-10) 

1 1^3n*^ * ‘^12''3n*^ ' ■ T5 * ^ o’*’ < ® o’ 

* !,2(»o”'3n«^’-'ll'W«o> 

(B-11) 

Solving Eq. (B-10) for Cj 2 yi®ids: 

r ~(i^k *''^’*'<^Po> - ‘=ll'*3n<^'’?o’ 

’*3o<"''So> 


(B-12) 
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Substituting Eq. (B- 12) into Eq. (B- 11) and solving for gives: 


=11 = I 


eo> 


1/ 




(B-13) 

where 


Plp< y = 'll<^>Y3n»»o) - 'l2<*'o>^3n«^> 

Once the parameters Y » v and n are specified and the second order 
functions are computed, all of the quantities in Eqs. (B-12) can then be 
calculated. Equation (B-12) is first solved for Cj 2 which is next substituted 
into Eq. (B-13) to obtain Cj|^. For v = 0.2, y = 1*^ and - ^ 23575 , 
the constants of integration become Cjj = 0.1346<1 and Cj 2 - “ 0.01078. 

After the constants and Cj 2 are obtained, Eq. (B-7) is used to 
calculate Fj(o), The functions and obtained by 

numerical integration. 

In computing the radial velocity by Eq. (43) the derivative f|(q' ) is 
also needed. Differentiating Eq. (B-7) yields; 

F '(o) 3 [ C j ^ - Ij2(a ) 1 a) + 3 [ Cj2 + Iii(o') ] Y3^(3v ) 


(B-14) 
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where and are simply the integrands in Eqs. (B>8) and (B<9). 


Solution for FjCw ) 

The third order function F^( «) is described by the inhomogeneous 
Bessel equation: 


5* k 




(B-15) 


where 




j^4( Y+1) ^ 2( - l)j R„f + 2( Y +ORnf' 


(B-16) 


The general solution of the homogeneous equation is 


’’5h<‘'> = Cjl V"> * Y„( 0.) 


(B-17) 


Adding the particular solution to Eq. (B-J7) yields the solution of the 
inhomogeneous equation as^^\ 


Fj(a) = [Cjj - a)l 3 n(a) + [ C52 + l5i( <^)1 Y^(o) 


(B-18) 
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where 


or 

f 3Jx)ft, 

vp„ 


(x)dx 


(R-I9) 


a 


152 (a) 


■/’ 

v^. 


Yj^(x) fe^(x)dx 


i^xr 


<ft-20) 


and W-( a) is the Wronskian given by: 


W3<^) = V«')Y^(a) -Y„(a)3^(a) 


(B-21) 


The constants of integration and C 52 must be determined 

using the boundary conditions given by Eqs. (35), Substituting Eq. (B^18) into 
Eqs. (35) yields: 


‘^5lV^o>*‘^52V^o>= [■^ 

*'52< Weo>-*5I< W9o> 


(B-22) 


(B-23) 


Unlike the corresponding equations obtained for f( a) and Fj(o ), 
Eqs. (B*22) and (B>23) contain an additional unknown parameter, the second 
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order eigenvalue correction P 2* additional relationship is needed 

to obtain ^ 2 * ^51 ^52* 

If ^2 known» the right hand sides of Eqs. (B-22) and (B-23) 

could be calculated. Denoting these quantities by r| and V 2 we have: 

V^'^o^ Sl^ V'^^o^S2 * ’’l 

(B-24) 

yPo>C5l" VPo>^2 = ^2 

This is a linear system of two equations in two unknowns which can be solved 
for C^i and C^2 determinant of the system is nonzero. The 

determinant A is given by 

^ = V'' - V '■ W ® o> ®-"> 

Comparing Eq. (B-25) with Eq, (8), which was used to determine shows 
that A =0, since is precisely the value that causes the right-hand- 
side of Eq» (B-23) to vanish. This implies that Eqs. (B-24) is either 
Inconsistent (i.e.» no solution exists) or that the system is indeterminate 
(i.e.y infinitely many solutions exist). The condition A ~ 0 also implies that 


(B-26) 


omm, r ^ 
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and the system is indeterminate if r^/r^ is also this ratio, otherwise no 
solution exists. This provides a method of determining ^2 requiring 
that solutions for and C ^2 exist, thus we require that: 

'•2 

or that 


('(vSo) - vej r;(vs„) 


[’f- - ej] - F5p(p„) 

■ W 


(B-27) 

where 

V Po> = 'Jl< » 0 >V? o> - >52< e o>V P o> <»■“> 

Solving Eq. (B-27) for ^2 introducing o = 3 ( v 0 ^) yields Eq. 

(36) which is used to calculate 62 * 

With ^2 given by Eq. (36), Eqs. (B-22) and (B-23) become an 
indeterminate system? that is, there are infinitely many solutions for C^j 
and C ^2 • 1 ^ order to obtain a unique solution for F^( a), an additional 
relationship between and C ^2 1*® found. This can be done if one 

states more precisely what is meant by the amplitude parameter e . 
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Until now» the amplitude parameter e was used in the perturbation 
analysis only as an ordering parameter and was not precisely related to a 
physical quantity such as pressure amplitude. In order to relate the 
theoretical calculations to experimental measurements, a more {.precise 
definition of e is needed. Since the Maslen and Moore perturbation scheme 
is based on the velocity potential, the amplitude parameter e is most 
conveniently related to the amplitude of the velocity potential. Accounting 
for the variation of frequency with amplitude, the expression for tp 
including all terms through third order is given by: 

cp s| eR^(a)+ e’^j^Fj(o)+ sin (t + n0 ) 

2 

+ - (y + lWa)] slu 2(t + ne) + c^F|(a)sin 3(t + ne ) 

(B-29) 

The parameter e cannot represent the zero-to-peak amplitude of cp since 

2 3 

this quantity includes terms proportional to € and e . For the same 
reason e cannot represent the amplitude of the first harmonic term since 

3 

this quantity includes terms proportional to e . A more restrictive 
definition of £ is therefore necessary. 

In the analysis to follow, e will be required to represent the 
amplitude of the term in the velocity potential with the time and space 
dependence of the first order solution, that is, Rj^(o) sin (t + n e). 
Considering the coefficient of sin(t + n0 ) in Eq. (B-29), there are three 
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terms: €R^(a), and The last term clearty does 

not have the same radial dependence as R^( (v) and does not need to be 
considered further. The term containing F^(<y) needs to be examined 
further since the homogeneous part of F^( or) is composed of the same 
linearly independent functions and Y^(«) that compose the radial 

acoustic eigenfunction R^(or ). This occurs because both ) and 

R^( or) are solutions of the same Bessel equation of order n. 


The homogeneous part of F^{ot ) can therefore be expressed as the 
sum of two parts, one of which is proportional to )• Thus ) can ^ 

be written as 

i Y_(vp ) . 

"sh*" > = j =5, - ^ j R„(a ) . Q„(« ) (B-30) 


where QAa) is a linear combination of 3 ) and Y (a), 
n n n 


Substituting Eq. (B-30) back into Eq. (B-29) yields the following coefficient 
of the acoustic solution R^(o)sin(t ne ): 


A = e + e 


^51" ar( vppS2 


(B-31) 


Since we desire that the amplitude of the acoustic solution A be equal to 

3 

the amplitude parameter e , the coefficient of e in Eq. (B-31) must 
vanish, which gives the following relation between and C^ 2 * 


'52 




'51 


(B-32) 
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This is the additionai relationship needed to uniquely determine F^(o). 
Equation (B-32) can be substituted into either Eq. (E-22) or Eq* <B-23) to 
eliminate one of the constants. For instance^ substituting Eq. (B-32) into Eq. 
(B-22) and solving for C^j yields 

Si = (B-33) 

For given values of Y 3 v , and n^ 8 2 calculated using Eq. 

(36). Then is computed using Eq. (B-33) and is substituted into Eq. 
(3-32) to determine C^ 2 * ^ Urst tangential spinning wavQ (n=l) and 

with Y = 1.<F, V = 0.2 one obtains ^2 “ ” •^9254, Cjj s ^ 0.0<>460 and 
Cj2= 0.01861. 

After the constants and C ^2 are determined* Eq. (B- 18) is used 
to caicuiate F^(o). The functions I^|(o) and determined by 

numerical integration. 

In computing the radial velocity by Eq. «^3) the derivative o) is 
also needed. Differentiating Eq. (B-18) yields: 

Fj(a ) = [ C^j - lj2^ 0 ^) ] «) + [ C^2 + *51^® ^ ^ 

- )3„(a ) + l5^(a )Y^(o ) (B-3<^> 

where ljj( ® ) and 1^2^ are simply the integrands in Eqs. (B-19) and 

(B-20). 
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Appendix C. 

Third Order Functions for Analysis with Steady Potentials 

Using the same techniques employed in Appendix the solutions of 
Eqs. (50) and (51) can be written as the sum of homogeneous and particular 
solutions as foliowst 






(C-1) 




(C-2) 


where 


I^l(a ) 




‘jr 


a 

_ f 

Jn(x)R^(x) 

11 
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”w(5y 


dx 


dx 
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Qjl(o') = 


Q^2(°^) s 


or 

f 

Y^(x)R^(x) 

4 

t{x) 

Oi 

f 


4 

X*W(«) 

or 

f 

Y„(x)R„(x) 


x^w(x) 


dx 


dx 


dx 


(C-6) 


(C-7) 


(C-8) 


and ) is given by Eq. (B-16) and the Wronskian W(a ) is given by 
W(a) = ya)Yj,(a)- (C-9) 


To complete the solution the constants of Integration C^|, C^ 2 » 


C ^2 determined. 


Applying the boundary conditions on F^(o ) given by Eqs. (57) yields: 


K 


‘^♦iVvPo) * V''^o> = ’*n<''Po’ <C-1I» 


SlVPo> * S2V»o> = 


mt 

^ *^n^ ^ o^ “ ^^9p^ ^ o^ 


(C-11) 
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where 


■= '«<WPo> 


(C-12) 


while the boundary conditions on F^(o ) given by Eqs. (Vf) yieldt 


J * C,,Y„(vSj = -2^ f'(v^ 


jlV'^o' ’ '-52'n"''’o' ■ ■ ''^2 




(C-iJ) 


<=5l’n<Po> ♦ C5 jY„(S„) = 


R'„< Po> 


“ ^o^ ♦ 2nKyFjy(p J (C-i^) 


where 


- '52<S o»n< »«> «:-»» 

•'5v< » ») = <55t< o> - V ^»n< «-'« 

Equations (C-13) and (C-14) contain four unknown quantities: the 
constants of integration and C^ 2 » second order eigenvalue 

correction ^ 2 , and the vortex constant K^. As shown in Appendix B the 


ORIGINU js 
OF POOS QUAtlTY, 


* ^ 4 


determinant of thi$ system is zero» therefore for solutions to exist the 
following relation must he satisfiedi 

■ VP J R'(vjl,) . 

*'0„) - h] ■*n<Po> ■ 


(Cvl7) 


Equation (C-17) reduces to Eq. (B-27) for ?: 0. If i$ known, 
Equation (C- 17) can be used to determine p2* method of 

Appendix B where c is required to represent the amplitude of the term in 
the velocity potential proportional to R^(cv)sin(t + n6), an additional 
relation between and C ^2 obtained. This relation is the same as Eq. 
(B-32), that is 

<=52 = (C.«) 

With p 2 Ky satisfying Eq. (C-17), Eq. (C-18) and either one of Eqs. 
(C>13) or (C^14) give three equations for the four unknown quantities 
Cj 2 » ^2 *^v* boundary conditions do not provide a method 

of uniquely determining ^2 p 5 (® ) when the vortex component of 

the acoustically driven steady flow is unknown. 
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Equations (C-10) and (C-U) are two relationships between the 
IS' Snowns C^^j, C^ 2 » *^s* determinant of this system is also zero, 

''.nd for solutions to exist the following condition must be satisfledt 








(C-19) 


In deriving Eq* (C-19), the steady potential source strength cancels 
leaving a relationship between fixed quantities which is probably not 
satisfied. Furthermore, even if Eq. (C-19) is satisfied, the boundary 
conditions do not provide a means for determining K^. 
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